Nodal quasiparticles and classical phase fluctuations in d-wave superconductors 
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We show that the nodal quasiparticles have significant effect on the classical phase fluctuations in 
a quasi-two-dimensional d-wave superconductor. They give rise to singularities in the temperature 
behavior of some of the coupling constants in the phase-only effective action. One of the consequences 
is that the classical XY-moAe\ is not adequate for the description of the superconducting fluctuations 
in d-wave superconductors at low temperatures. 
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In superconductors with low superfluid density ps , the 
' phase fluctuations of the order parameter can play a sig- 
nificant role. Emery and Kivelson I'] suggested that it is 
the classical phase fluctuations of the pre-formed Cooper 
pairs that are responsible for the pseudogap phenomenon 

■ in the underdoped high-T^ materials (HTSC's) [3|. This 
[ idea has been further elaborated by many authors, see, 

e.g., the review article . On the other hand, it was ar- 
I gued in Refs. that the phase fluctuations in HTSC's 
are important down to T <C Tc, which could explain a 
linear in T behavior of the magnetic penetration depth 
A(T) at low temperatures |^, even for a nodeless order 
parameter (in a superconductor with lines of nodes, a lin- 
ear T-dependence of A occurs already at the mean-field 
level due to the thermally excited nodal quasiparticles 

: 0). 

The effects of the phase fluctuations in HTSC's are 
usually studied using either the classical or the quantum 
versions of the XY-model, which is believed to provide a 
correct description of any system with a one-component 
complex order parameter of fixed magnitude. The energy 
of the classical two-dimensional (2D) XY-mode\ in the 
[ absence of external fields has the form 

■ £jfY =^ Jrr'[l-COs(0r-^r')] (1) 

rr' 

] Here Jrr' are the phase stiffness coefficients, and 9r is the 
phase of the order parameter at site r of a coarse-grained 
square lattice. The lattice spacing d is of the order of 
the superconducting correlation length ^, which is almost 
constant at low temperatures. The XY-mode\ has also 
been applied to other systems, such as fabricated arrays 
of Josephson junctions, or granular superconductors (see, 

I e.g., Ref. and the references therein), where d is equal 
to the distance between grains. 

■ Although the simplicity of the model ^ is physically 
appealing, it should be stressed that its rigorous micro- 
scopic derivation for homogeneous high- Tc superconduc- 
tors does not exist. The usual way of justification (see, 
e.g., Ref. 0) involves expanding the cosine in the con- 
tinuous version of Eq. and matching the coefficients 
with those in the Gaussian phase-only action obtained by 
integrating out the fermionic degrees of freedom. At low 
temperatures, varies slowly at the length scale of d and 
vortices can be neglected. Then, assuming the nearest- 



neighbor interaction Jrr' = J in Eq. , introducing the 
superfiuid velocity Vs = V6'/2m (m is the electron mass, 
and we use such units that h — c = e — 1) and expanding 
the cosine in powers of and its gradients, we obtain 

+ f «.+<,) + ■■■}■ (2) 

Here ps — Am?, J is the mass density of the superfluid elec- 
trons, K = rn^d^J, and B = -{8/3)m'^d'^J. The ellipsis 
stands for the terms of higher order in Vs or its gradients, 
such as (Vvs)'^V3, etc. Given J and d in Eq. one can 
easily work out the coefficients in all orders of the expan- 
sion. [If to take into account the next-nearest-neighbor 
interactions Jrr' 0|, then the tensor structure of the 
gradient and interaction terms in Eq. will change, 
for example the terms proportional to xVs,x){'^yVs,y), 
or v'^Vy will appear.] An important feature of the XY- 
model (0) is that K, B and all the other expansion 
coefficients are temperature-independent at low T . 

That the consistent microscopic theory fails to repro- 
duce the quantum generalization of the XY-raodel has 
already been noticed in Ref. this Letter, we show 

that the long-wavelength field theory even for the classi- 
cal phase fluctuations in a c?-wave superconductor at low 
temperatures does not have the form of Eq. ((2Jl . Briefly, 
the microscopic theory shows that the coefficient in front 
of the gradient term in the energy of fluctuations diverges 
at r — !■ 0, indicating a strongly non-local interaction of 
the phase fluctuations. In addition, the Taylor expansion 
of the energy of the phase fluctuations in powers of the 
uniform superfluid velocity breaks down at T ^ 0, lead- 
ing to a non- analytical dependence on Vg. To the best 
of our knowledge, the non-local and non-linear effects on 
the phase fluctuations in high-Tc superconductors have 
not been emphasized before. 

The starting point of our microscopic analysis is a 
tight-binding Hamiltonian on a 2D square crystal lattice 
with the "microscopic" lattice constant a: 

^ = E^ij^U'^-5E4^y- (3) 

ij (ij> 

Here = — — /if^y, iij is the hopping amplitude, p 
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is the chemical potential, g > is the superconduct- 
ing coupling constant, and the operator iJy = (cj|Ci| — 
CjiCii)/V2 creates a singlet pair of electrons at the 
nearest- neighbor sites (ij). The sites of the crystal lattice 
are labelled by i, in contrast to the sites of the coarse- 
grained lattice in Eq. (QJ . We assume zero external mag- 
netic field and neglect disorder and the electron-electron 
interactions other than that responsible for the Cooper 
pairing (it can be shown that the results below for the 
classical phase fluctuations are not affected by the latter 
assumption) . 

Following the standard procedure, the partition func- 
tion for the model Q can be written as a Grassmann 
functional integral. The interaction term are then de- 
coupled using the Hubbard-Stratonovich transformation, 
and the fermionic fields are integrated out. After these 
manipulations [T^ . we end up with a representation of 
the partition function as a functional integral over a com- 
plex bosonic field A which is non-zero only on the bonds 
between the nearest neighbors: 

Z = J PAij (r)I? A?j (r ) 6"^=^^ ^"^"'^l (4) 

Here < r < /3 is the Matsubara time, P = 1/T, and the 
effective action is given by 

5,// = -Trln(?-i + - /"dr ^|Ay|2, (5) 

9 Jo 

with 

We use the notation Tr for the full operator trace with 
respect to the space-time coordinates and the Nambu 
indices, reserving the notation tr for the 2x2 matrix 
trace in the electron-hole space. 

The mean-field BCS theory corresponds to a station- 
ary and uniform saddle point of the functional integral 
Q. The d-wave saddle point in the momentum rep- 
resentation is given by Ak = Ao(r)0k, where (/)k = 
2(cos kxCL — cos kyO) is the symmetry factor. The temper- 
ature dependence of the order parameter is determined 
by the standard BCS self-consistency equation, general- 
ized for the case of an anisotropic order parameter (see, 
e.g., Ref. 01). Deviations from the mean-field solution 
can be represented in terms of the amplitude and phase 
fluctuations: Aij(r) = [A?^ -f (5Aij(r)]e^'^'j(^), where aI"' 
is the Fourier transform of Ak. We neglect the amplitude 
fluctuations because they are gapped and therefore make 
a negligible contribution at low temperatures. Since the 
number of bonds in a square lattice is twice the number 
of sites, one should introduce two phase fields on sites, 
6*1 (t) and 0i(r), to describe the phase degrees of freedom. 
One possible parametrization is = 6i for rj = ri + ai, 
and = Oi + Oi for rj = Vi -\- ay. As shown in Ref. 
, the fiuctuations of 0, which describe a change in the 



symmetry of the order parameter from d-wave to d -|- is- 
wave, can be neglected. Making a gauge transformation 
^ij' ^ ^ij' = UlG^^U^, where U = exp[-(i/2)0T3] (tj 
is the Pauli matrix in the electron-hole space) and 
expanding the Tr In^"^ = Tr \nQ~^ in powers of the 
small temporal and spatial gradients oi 9, we obtain a 
phase-only effective action. 

Motivated by the arguments of Emery and Kivelson 
that the phase fluctuations in HTSC's are predominantly 
classical down to low temperatures, we consider in this 
Letter only the classical case, which corresponds to ne- 
glecting the r-dependence of 0. The effective action then 
becomes Seff = Pi£o + £), where Sq is the mean-field 
energy, and £ is the energy of phase fluctuations. In the 
Gaussian approximation, 

^ = \ j d^'>'id^r2 ICijiri - r2)t;s,j(ri)t;sj(r2), (7) 

where i, j = x,y, and the Fourier transform of the kernel 
/C(R) has the form 

/C.y(q) = PsfiSij - m^ny(q). (8) 

To simplify the notations, from this point on we assume 
a quadratic band dispersion ^k = (fc^ — kp)/2m* with 
the effective mass m* equal to the electron mass m. In 
Eq. (O, psfi — nom/a^ is the superfluid density at zero 
temperature, no is the average number of electrons per 
site, and 11^ is the current-current correlator: 

^ f d^k 
n,,(q) = -2T^ / 7^-^^^- 

xtr[go(k + q,cj„)e^o(k,w„)]. (9) 

Here t/o(k,a;„) = {iutnTQ — ^i^t^ — AkTi)^^ is the mean- 
field matrix Green's function, ujn — (2?! -I- l)7rT is 
the fermionic Matsubara frequency, and v = Vk^k is 
th e Fermi velocity. The energy of excitations, i?k = 
•\/^^ + AJ, vanishes at the four gap nodes k = k„ = 
kF{±x±y)/V2, n = lH-4. 

At low temperatures, only the nodal excitations with 
k = ki -|- (5k are important, which allows one to linearize 
both the quasiparticle dispersion ^k — vpSk^ and the 
order parameter Ak — fA<5fc||, where Sk± and are 
the components of i5k perpendicular and parallel to the 
Fermi surface. The anisotropy ratio vf/va depends on 
the material and the doping level, e.g. vp/vA = 14 and 
19 in optimally doped YBCO and Bi-2212, respectively 
|15| |. The fermionic excitations have a conical spectrum: 
Sk = y^w|(5fc2 + u|(5fcj^ HE IS- Calculating the Mat- 
subara sum in Eqs. (jSJ), we find 

ICxxiq.) = K^vvifi) = Ps.o - ^+(q) 

(10) 

where 
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Here 



7i.2(q) = ^ 



?{qx ± Qy)^ + vliq^ T Qy)'^ (12) 



are the energies of the nodal quasiparticles with (5k ~ q, 
and s{x) is a dimensionless scaling function, which is de- 
fined by a rather cumbersome integral, whose asymptotic 
expansions at large and small x are 



six) 



21n2 + a;V24 
Trx/8 



at X 
at X 





oo. 



(13) 



We would like to emphasize that the expressions (|ll|l are 
exact in the nodal approximation, i.e. for the conical 
quasiparticle spectrum. In terms of q, the applicability 
region of the nodal approximation is 71,2 (q) <C Aq. 

At finite temperatures, the long-wavelength behavior 
of the kernel /C(R) corresponds to the limit 7i.2(q) ^ T. 
In this case, according to Eqs. (|10|l and H13|) . the Fourier 
components /Cij(q) are analytical functions of momen- 
tum, which can be written in the following form: 



K.xx{ci) = K^yy (qi) 



(14) 



)Cxy{(l) = ICyxiq) = -2psl'^qxqy, 



where Ps{T) — ps.o — {2\n2 j 'r:){m?vp /v/^)T is the mean- 
field superfiuid mass density suppressed by the thermally 
excited nodal quasiparticles [l3- To leading order in 
(va/vf)^ <C 1, the characteristic length ^ is given by 



^(T) = ^— 

\A8ttvaPsT 



1/2 



C 



1/2 



Co, 



(15) 



where Tc is the critical temperature, ^0 = vf/2ttTc is the 
BCS coherence length, and C = [Trm^TcVF I^2psVaY/'^ 
is almost constant at low T. For a circular Fermi surface 
and T = 0, C = ^ [t:^ /12){vf /va){Tc/ojf) ~ 1. 

The quadratic q-dependence (I14|l implies that the ker- 
nel /C(R) in real space is proportional to exp(— |R|/^). 
This behavior is similar to that of the electromagnetic 
response function in conventional s-wave superconduc- 
tors, see e.g. Ref. 18] . An important difference however 
is that the characteristic length of fluctuations ^ is now 
temperature- dependent: ^(T) ~ T~^l'^ at T — > 0. Quali- 
tatively, this can be attributed to the partial filling of the 
gap nodes by the thermally-excited quasiparticles, which 
creates an effective temperature-dependent s-wave gap. 
Note that the length ^, defined by the gradient expan- 
sion of the current-current correlation function, Eq. (|14|l . 
is different from other characteristic lengths discussed in 
the literature: the coherence length ^0, which is the cor- 
relation length of the amplitude fiuctuations, and the size 
of a Cooper pair ^pair^ which is infinite in the d-wave case 
|l9j| . In a conventional s-wave superconductor, all three 
lengths are of the same order. 



If the typical length scale of the phase fluctuations ex- 
ceeds C(r), then the energy of such fluctuations at finite 
temperatures in the Gaussian approximation can be writ- 
ten asf = (l/2)(/ d2i?/C,,)(/ d^r v^) = (p,/2) / v^, 
which indeed reproduces the quadratic term in the con- 
tinuum version of the A"!"- model (|2}. However, compar- 
ing the gradient terms in the A"F-model with the 
microscopic result (|14f) . we see that they are clearly dif- 
ferent. First, the coefficient K turns out to be singular 
at T ^ 0. Second, the sign and the tensor structure of 
the gradient terms are different. 

Because of the divergence of ^, the gradient expansion 
of £ breaks down at T ^ 0. As evident from Eqs. H10() 
and H13|) . at T = the matrix elements /Cy (q) are non- 
analytical functions of momentum. One can therefore 
expect that the kernel /C(R) decays algebraically at large 
R. Indeed, calculating the Fourier integrals, we obtain 



/C±(R) 



IGttw^ 



[/+(R)±/_(R)], (16) 



where 

/±(R) = 



Rx i Rv\^ f RxT R-i^ 



^ x J-i-y 
Vp 



-1 -3/2 



(17) 



Thus, the functions /Cy (R) are strongly anisotropic and 
decay slowly, as |R|^^, along a given direction. These 
expressions also describe the behavior of ICij (R) at finite 
T and at shorter distances, such that T ^ 7i,2(q) *C Aq. 
At larger distances there is a crossover to the exponential 
decay obtained above. 

In order to go beyond the Gaussian approximation, 
one has to take into account the higher-order terms in 
the Taylor expansion of the effective action |SJ) in powers 
of V0 or Vg. Assuming a uniform superfiuid velocity, one 
can derive from (jSJ the following exact expression for the 
energy density of fluctuations valid in all orders in 
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V 



= -T 



(2^ 



cosh -^jr- 

In 221 



cosh ■ 



2T 



cosh 



2 Bk 
2T 



(18) 



where V is the system volume, ii^k,± = (^k,+ — ^k,-)/2± 
\/ i^k,+ + ^k,~)^/4 -I- are the excitation energies af- 
fected by the superflow via the Doppler shift in the en- 
ergy spectrum, ^k,± — fkimv,- Keeping the first terms 
in the Taylor expansion of £, we have 



£_ 

V 



Ps 



4 ^ " 



4w^ ) 

s,x s.y } 



(19) 



where Ps(T) is the mean-field superfiuid density defined 
above, and B{T) = const - 2m?(^{T)ps oc T^^, to lead- 
ing order in T~^, with f given by Eq. (|15|1 . We see that 
this expression does not have the form of the continuum 
Xy-model since the coupling constant B is divergent 
at r ^ 0. The singularity of B in Eq. (|19|) indicates that 
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the Taylor expansion of the energy of fluctuations £ in 
powers of breaks down at low temperatures. Indeed, 
at T = the k-integral in Eq. H18|l can be calculated, 
with the result 



3, ,2 



m V 



_ ^ Ps,0 2 ^_ 

y ~ 2 12V27rwA 



{\Vs 



|3 



(20) 



It is this non-analytical dependence that is responsi- 
ble for the non-linear Meissner effect in HTSC's 
[to see this, one should use the gauge-invariant super- 
fluid velocity in the presence of a vector potential A, 
V, ^ = {l/2m)Ve - A]. 

The non-linear terms in £ are responsible for the in- 
teraction of the thermally-excited classical phase fluctu- 
ations, which in turn will lead to a renormalization of the 
bare superfluid density ps (T) . Despite the divergence of 
the uniform coupling constant _B at T ^ 0, it is possi- 
ble to calculate the lowest-order correction to ps, which 
turns out to be non-singular (and moreover it identically 
vanishes at T = if the quantum effects are taken into 
account 113). The formal reason is that the contribution 
of the uniform vertex vanishes because of the momentum 
structure of the corresponding diagrams. 

In conclusion, we have shown that the effect of the d- 
wave nodal quasiparticles on classical phase fluctuations 
cannot be ignored. The microscopic derivation of a low- 



temperature effective energy of the classical phase fluctu- 
ations gives a result that is very different from the XY- 
model. The nodal quasiparticles with conical dispersion 
lead to the appearance of a new, temperature-dependent, 
length scale ^, whose divergence at T ^ causes both 
the gradient expansion and the Taylor expansion in pow- 
ers of Vg of the energy of fluctuations to break down. 
The physical origin of the first effect is the same as that 
of the non-local Meissner effect in d-wave superconduc- 
tors since the superconducting order parameter Ak 
has nodes on the Fermi surface, then close to the nodes 
anisotropic coherence length Wi?/|Ak| exceeds the Lon- 
don penetration depth Aq, and the local electrodynam- 
ics breaks down. The failure of the Taylor expansion is 
closely related to the non-linear Meissner effect l2^: the 
superfluid velocity acts as a pair-breaker in d-wave super- 
conductors, creating a finite density of normal excitations 
even at T = 0. The backflow of these excitations leads to 
a non-analytical dependence of the electromagnetic re- 
sponse functions on v^. The divergence of ^ at T — + 
might be cut off by the effects not included in our analy- 
sis, such as quantum fluctuations, inter-layer tunnelling, 
or disorder. 
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